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Abstract. With a collecting area of 70 000 m^, the Five hundred metre Aperture Spherical Telescope (FAST) will 
allow for great advances in pulsar astronomy. We have performed simulations to estimate the number of previously 
unknown pulsars FAST will find with its 19-beam or possibly 100-beam receivers for different survey strategies. 
With the 19-beam receiver, a total of 5 200 previously unknown pulsars could be discovered in the Galactic plane, 
including about 460 millisecond pulsars (MSPs). Such a survey would take just over 200 days with eight hours 
survey time per day. We also estimate that, with about 80 six-hour days, a survey of M31 and M33 could yield 
50-100 extra-Galactic pulsars. A 19-beam receiver would produce just under 500 MB of data per second and 
requires about 9 tera-ops to perform the major part of a real time analysis. We also simulate the logistics of 
high-precision timing of MSPs with FAST. Timing of the 50 brightest MSPs to a signal-to-noise of 500 would 
take about 24 hours per epoch. 
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1. Introduction 

FAST, the Five-hundred-metre Aperture Spherical 
Telescope, is an Arecibo-style telescope currently under 
construction in China. FAST will be located in a karst de- 
pression in Guizhou province at a latitude of about +26° 
and i s planned to be completed around 2014 (|Nan I2OO6L 



200i). With an illuminated aperture of 300 m diameter, it 



will be an order of magnitude more sensitive than the 
steerable 100-metre telescopes at Effelsberg and Green 
Bank and about twice as sensitive as Arecibo. Moreover, 
FAST will cover a frequency range of 70 MHz to 3 GHz and 
it will be able to observe at a zen ith angle of 40 degrees 
without a significant loss in gain ( Jin et al. 20081 ). By il- 
luminating only part of the 300-metre aperture the zenith 
angle can be extended even further; at a zenith angle of 60 
degrees about half the aperture is still illuminated. With 
its sensitivity and its ability to see the majority of the 
sky, including a large part of the Galactic plane, FAST 
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has excellent po tential for finding radio pulsars (see also 
iNan et al.|[2006h . 

Radio pulsars provide insights into a rich variety of 
physics and astrophysics . Apphcations to date (see e.g. 
iLorimer fc Krarrier 20051 ) include the study of the Milky 
Way, globular clusters, the evolution and collapse of mas- 
sive stars, the formation and evolution of binary systems, 
the properties of super-dense matter, extreme plasma 
physics, tests of theories of gravity and as cosmologi- 
cal probes. A goal for the coming decade is the direct 
detection of low-frequency gravitational waves via high- 
precision timing of an array of millisecond pulsars (MSPs) 
( Jenet et al ] l2005h . Such an array would be greatly fa- 
cilitated by further discoveries and timing of MSPs with 
FAST. 

The aim of this paper is to investigate the possibili- 
ties of FAST for finding and timing radio pulsars. In §2 
we present a simulation of a FAST survey of the pulsar 
sky, where we look into the survey speed and the num- 
ber of pulsars that can be detected for different receivers. 
In §3 we present calculations to estimate the data rates 
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Table 1. Expected system parameters of FAST and the 
pulsar survey. 



Parameter 


Value 


Gain 


16.5 KJy-' 


System temperature 19-beam receiver 


20 K 


System temperature 100-beam receiver 


30 K 


Centre frequency 


1315 MHz 


Bandwidth 


400 MHz 


Number of frequency channels 


9 500 


No. of polarisations 


2 


Beam FWHM 


3'.4 


No. of beams 


19 or 100 



and computational requirements for a pulsar survey with 
FAST. In §4 we look into the possibihty for FAST to par- 
ticipate in a global pulsar timing array effort. §5 contains 
a discussion of the results, as well as a comparison between 
FAST and the Square Kilometre Array. 



2. Survey simulation 

In the past, the southern hemis phere has been the sub- 
ject of extensive pulsar surveys (IManchester et al. 119961: 



Lvne et all [l998l: 



Manchester et al.M200lD . FAST will be 



able to complement these surveys in the northern hemi- 
sphere with an even greater sensitivity. Initially, FAST will 
have a 19-beam receiver with a frequency range of 1.1- 
1.5 GHz and a system temperature of 20 K (not including 
the sky temperature). A possible future p hased a rray feed 
(FAF) is planned to have over 100 beams (jjin et al.,.2008 ) 
with the same frequency range and a system temperature 
of 30 K. Table [T] shows the parameters involved in the sur- 
vey simulation. The number of beams will determine the 
speed of the survey and the dwell time per pointing. We 
therefore perform the simulations as a function of both 
frequency and field of view (FoV) . 

2.1. Simulation method 

We performe d Mo nte Carlo simulations following 
Lorimer et al. I (l2006h using the psrpojj^ packa g e (cf . 



Smits et al.ll2009t ). In their study. iLorimer et al.l (|2006f) 



used the results from recent surveys with the Parkes 
Multibeam system to derive an underlying population of 
pulsars with an optimal set of probability density func- 
tions for pulsar period (P), 1400-MIIz radio luminos- 
ity (L), Galactocentric radius (i?) and height above the 
Galactic plane (z). We make use of these results in our 
simulations de scribed below which u se as a starting point 
model C' from 'Lori mer et aL ( 20061 ). 

Our simulation procedure begins by generating a pop- 
ulation of normal pulsars which beam towards the Earth. 
Each pulsar is assigned a value of P, L, R and z based on 
the assumed probability density functions. For the distri- 
butions in P, R and z, we use the distributions from C' of 



Lorimer et al. ( 20061 ). To com pute L, we adopt the log- 
norma l distribution found by iFaucher-Giguere fc Kaspi 
(|2006l) . This has the advantage of not requiring a spe- 
cific lower bound in L, as is th e case for the power- law 



luminosity models considered bv lLorimer et al. I (|2006h . 

We compute intrinsic pulse widths using the follow- 
ing self-consistent approach. Assuming a simple geome- 
try with circular beams of radius p, the pulse width can 
be found from the inclination angle between the pulsar 
spin and magnetic axis, a, and the impact parameter be- 
tween the magnetic axis and the line of sight, /3. Following 
Kramer et al. (1998), we use the empirical relationship be- 
tween p and spin period P (s) defined as follows: 



Pmodcl 



5 4op-i/2 if p > 30 ms 
31.2° otherwise, 



where the angle of 31.2° is the value of the above ex- 
pression evaluated for P = 30 ms. To model the scatter 
about this relationship (see, e.g.. Fig. 3.5 from Lorimer & 
Kramer 2005), we use a Monte Carlo approach to find a 
dithered beam radius 



P 



(1) 



where p is a random number drawn from a flat distribution 
in the range — 0.15 <p< 0.15. This, somewhat arbitrary, 
scheme dithers the beam radius logarithmically by a third 
of a decade and provides a good match to the observed 
scatter. The pulsar is assumed to be beaming towards us 
and the impact parameter (3 is chosen from a simple flat 
distribution in the range —p < P < p. Finally, we as- 
sume that the magnetic inclination angle a is randomly 
distributed, that is 



a = arccosq, 



(2) 



where g is a random number drawn from a flat distribution 
in the range Q < q < 1. With these quantities defined, 
the observed pulse width W is found using the following 
geometrical relationship (see, e.g., Gil et al. 1984): 



sm 



sm 



sm 



(f) 



sin a • sin(a + (3) 



(3) 



^ |http: //psrpop. sour ceforge .net] 



To compute the expected DM and scatter broadening 
effects on each pulse, we use th e NE2001 electron den- 
sity model (jCordes fc Lazioll2002l ). Note that, since we are 
primarily concerned with the distant population of highly 
dispersed pulsars in these simulations, we do not attempt 
to account for interstellar scintillation. Finally, to allow 
us to extrapolate the 1400-MIIz luminosities to other sur- 
vey frequencies in the next step, we make the reasonable 
assumption t hat pulsar spectra c an be approximated as 
a power law ( Lorimer et al. 19951 ) and assign each pulsar 
a spectral index drawn from a normal distribution with 
mean of -1.6 and st andard deviati o n 0.35 . 

As described by ILorimer et aL ( 20061 ). when generat- 
ing the normal pulsar population according to the above 
criteria, we keep track of the number of pulsars detectable 
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by the detailed model of the Parkes Multibeam Pulsar 



10000 



Surveys of the Galact ic plane ([Manchester et alj l200lh 
and at high latitudes ( Burgav et alJ 2006^ . Oursimula- 
tions terminate when the numbers of model detectable 
normal pulsars matches the 1005 pulsars found by these 
surveys. This results in an underlying sample of 120,000 
model pulsars. 

To model the MSP population, we follow the same 
approach and assumptions as for the normal pulsars de- 
scribed above, but with two exceptions: (i) model z heights 
are chosen from an exponential distr ibution with a mean 
of 500 pc (' Cordes fc Chernofilll997[) : (ii) the underlying 



period distribution used is taken from a recent study of 
the MSP population (Lorimer et al. 2009; in preparation). 
Note that implicit in this approach is the notion that 
the luminosity functions for normal pulsars an d MSPs are 



identi cal; this is consistent with earlier results iLvne et al 
(Il998h . 



In addition to the Parkes Multibeam Surveys con- 
sidered above for the normal pulsars, we also make use 
of the Parkes 70-cm pulsar survey of th e southern sky 
( Manchester et al. 19961 : Lvne et al. 1998 ) and two inter- 
mediate latitu de surveys carried out by the Parkes mult i- 



beam system (jEdwards et al.l l200lt iJacobv et al.l 120071) 



Our simulations are normalized such that they terminate 
when the total of 49 MSPs detected by these surveys is 
reached. This results in an underlying sample of 23 000 
model MSPs. 

Once the pulsar population was determined, we per- 
formed simulations to find the number of pulsars FAST 
would detect as a function of frequency. For this purpose, 
we allowed the centre frequency to range from 400 MHz 
to 1.4 GHz. The bandwidth was kept at one third of the 
centre frequency, the observation time per pointing was 
kept constant at 600 s and the gain was kept constant 
at 16.5 KJy~^. The frequency-dependent sky temperature 
was added to the system temperature. Considering the 
beamwidth of about 3.4 arcminutes at 1.4 GHz, we lim- 
ited the surveys to Galactic latitudes of either between 
±5° or between ±10°. Further, we simulated the number 
of pulsars FAST will detect as a function of observation 
time. All detected pulsars include already known pulsars. 

2.2. Simulation Results 

Fig. [1] shows the number of normal pulsars and MSPs that 
were detected in the simulation of the Galactic plane as 
a function of centre frequency. The figure shows that for 
both normal pulsars and MSPs, the number of detected 
pulsars is reasonably stable from centre frequencies of 1 
to 1.4 GHz. Similar results are found for observation times 
of 120 and 1 800 seconds. The number of detected pulsars 
as a function of observation time is shown in Fig. [2j This 
figure shows that, while the number of detected pulsars 
continues to climb with increasing observation time, the 
rate of increase slows significantly after times of about 600 
s. 
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Normal pulsars with |b|<5 
Normal pulsars with |b|<10 
IVIiliisecond pulsars with |b|<5 
Millisecond pulsars with |b|<10 
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Fig. 1. Number of normal pulsars and MSPs detected in 
simulations of FAST surveys of the Galactic plane as a 
function of observation frequency. The surveys were lim- 
ited in Galactic longitude by 20° < 1 < 90°. The sohd hues 
and the dotted lines show the number of normal pulsars 
and millisecond pulsars, respectively. The squares indicate 
a survey region limited in Galactic latitude by |b| < 5°. 
The circles indicate a survey region limited in Galactic 
latitude by |b| < 10°. The observation time was kept con- 
stant at 600 s and the bandwidth was kept at one third of 
the observation frequency. 
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Fig. 2. Number of normal pulsars and MSPs detected in 
simulations of FAST surveys as a function of observation 
time per pointing. The surveys were limited in Galactic 
latitude and longitude by |b| < 10° and 20° < 1 < 90°. 
The solid line and the dotted line show the number of nor- 
mal pulsars and MSPs, respectively. The centre frequency 
is 1.315 GHz with a bandwidth of 400 MHz. 



The 19 beams in the multibeam receiver of FAST have 
a total FoV of about 0.061 deg^. A possible future PAF 
with over 100 beams will have a total FoV of 0.32 deg^. 
Table shows the total survey time for different obser- 
vation times per pointing, for different survey regions for 
both a 19-beam and a 100-beam receiver. The total survey 
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Table 2. Total survey time for different FAST surveys. Except for model 7, the survey region is limited in Galactic 
longitude by 20° < 1 < 90°. The survey time is expressed in number of days with eight hours of observation. The 
number of detected normal and millisecond pulsars are rounded to the neares t 100, and 10, res pectively. The number 
of known pulsars has been obtained from the ATNF Pulsar Catalogue (Manc hester et al. 2005f ). 



Parameter 


Model 1 


Model 2 


Model 3 


Model 4 


Model 5 


Model 6 


Model 7 


Survey region 


|b|<5° 


|b|<10° 


|b|<5° 


|b|<10° 


|b|<5° 


|b|<10° 


Visible Sky 


Pointing time (s) 


600 


600 


600 


600 


1800 


1800 


120 


Number of beams 


19 


19 


100 


100 


100 


100 


100 


Number of normal pulsars detected 


5100 


6 300 


4 200 


5 200 


5 700 


7000 


4 500 


Number of MSPs detected 


470 


670 


380 


550 


550 


770 


500 


Currently known normal pulsars in region 


352 


418 


352 


418 


352 


418 


662 


Currently known MSPs in region 


14 


20 


14 


20 


14 


20 


62 


Total survey time (eight hour days) 


229 


454 


44 


86 


131 


259 


296 



time is given in days assuming eight hours of observation 
time per day. With the initial 19-beam receiver about 5 200 
previously unknown pulsars, including 460 MSPs, can be 
found by searching the region with 20° < 1 < 90° and 
|b| < 5°, using 10-minute pointings. This survey would 
take just over 200 days to complete. Using the 100-beam 
receiver to survey the same region with 30-minute point- 
ings, 5 900 previously unknown pulsars can be found in 
130 days, including about 540 MSPs. Doubling the survey 
area to |b| < 10° will add about 1 100 pulsars using 10- 
minute pointings and about 1 450 pulsars using 30-minute 
pointing. The 100-beam receiver will also enable an all- 
sky survey. Using 2-minute pointings, 4 300 previously un- 
known pulsars can be found in 300 days, including about 
440 MSP's. 

Fig. [3] shows the distribution of detected pulsars for 
all the models. The top two plots show the normal pul- 
sars that are detected in the models that are limited to 
the Galactic plane. The plots directly below show the mil- 
lisecond pulsars that are detected in the same models. The 
bottom plot shows the all-sky survey. 

2.3. Survey of M31 and MSB 

The great sensitivity of FAST will permit searches for 
pulsars in other galaxies. We can estimate the number 
of pulsars that FAST can detect in a survey of M31 and 
M33, both visible for about six hours per day from FAST. 
Table [3] shows the parameters of these galaxies. The min- 
imum detectable flux density of a pulsar in a survey is 
given by: 



= 20K, Tsky = 1.4K for M31 and Tsky = 1.1 K for M33, 



S/N (Tsy, + Tsky) 
Gy^UptobsB 



W 



P-W 



mJy, 



(4) 



where S/N is the signal-to-noise ratio required for a de- 
tection, Tsys is the system noise temperature (K), Tgky is 
the sky temperature (K), G is the gain of the telescope 
(KJy~^), rip is the number of polarisations, iobs is the 
observation length (s) per pointing, B is the observation 
bandwidth (MHz), P is the period of the pulsar (s) an d W 
is the equivalent width (s) ( Lorimer fc Kraniieill2005 ). For 
an observation frequency of 1.315 GHz and S/N = 9, Tgys 



G = 16.5 KJy- 



= 2, tobs = 6 hours, B = 400 MHz 



and assuming an equivalent width of 5% of the period 
we obtain a minimum detectable flux density of 600 nJy. 
This leads to a minimum detectable luminosity of 360 and 
560mJykpc^ for a pulsar in M31 and M33, respectively. 
There are at least 55 pulsars known in our Galaxy with 
a higher luminosity at 1.4 GHz than 360mJykpc^ and 32 
pulsars higher than 560mJykpc^. If we correct for the 
mass ratio between our Galaxy and the galaxies M31 and 
M33, then we can expect of order 80 pulsars to be de- 
tected by a FAST survey of M31 and about 3 pulsars in 
M33. However, the supernova rate per unit mass in our 
Galaxy is thought to be higher than in M31 and consider 



ably lower than in M33 (jBerkhuii seni 1 1 984t iGordon et al 
1998). Also, since the angular extent of M33 is about 1 
quarter that of M31, the number of pulsars found per 
pointing might not differ much between M31 and M33. We 
estimate that the total number of pulsars found in a FAST 
survey of both M31 and M33 would lie between 50 and 100. 
Finding such a number of pulsars would allow studies of 
the pulsar population in these galaxies and would provide 
a valuable probe of the inter-galactic medium. The com- 
bined angular extent of M31 and M33 is 3.8 deg^, which 
means about 1 500 beams are required. With a 19-beam 
receiver, this would require 79 pointings each of duration 
six hours, leading to a total of 470 hours of observing time. 

3. Data processing 

The total amount of data from a FAST pulsar survey is 
given by 



total 



j-y, ^survey i» bits 
-Dinh.^ iV„, 



FoVbc 



Vpo/- 



Bytes, 



(5) 



where B is the bandwidth, tots is the observation time 
per pointing, FoVbeam is the FoV of 1 beam, fisurvcy is the 
total field of the survey, iVpoi is the number of polarisations 
and iVbits is the number of bits used in the digitisation. 
Figure |4] shows the total amount of data from a FAST 
pulsar survey as a function of the total survey field for 
different observation times. We have assumed a bandwidth 
of 400 MHz, a FoVbeam of 0.0025 deg^, 1 polarisation and 
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PSR's in 20" < 1 < 90°, Ibl < 5° using 600 and 1800 s pointing time 




1 800 s pointing time 
600 s pointing time 
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PSR's 20° < 1 < 90°, Ibl < 10° using 600 and 1800 s pointing time 




1800 s pointing time 
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MSP's 20 < I < 90 , Ibl < 10 using 600 and 1800 s pointing time 
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600 s pointing time 
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Fig. 3. Distribution of detected normal pulsars and MSPs for different models. The top left plot shows the distribution 
of normal pulsars from models 1 (similar to model 3) and 5. The top right plot shows the distribution of normal pulsars 
from models 2 (similar to model 4) and 6. The two plots in the the middle show the distribution of millisecond pulsars 
from the same models as the plot directly above. The bottom plot shows a Hammer- Aitoff projection of the distribution 
of normal pulsars and millisecond pulsars from model 7 in galactic coordinates. The black and white striped boxes 
mark the areas defined by 20° < 1 < 90°, |b|<5° and |b|<10°. 



2 bits per sample. Table [4] shows the total amount of data 
for the seven different survey models from table [2] 



The data acquisition rate from a pulsar survey sam- 
pling the detected outputs of a filterbank of total band- 
width B and channel bandwidth Av is given by: 

'Drate = A^boams ^ BpS, (6) 

Esamo ° 
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Table 3. Parameters of M31 and M33. 



Galaxy 


Declination 


Distance (kpc) 


Mass (M0) 


Angular extent 


M31 
M33 


41° 16' 
30° 39' 


772±44 fRibas et al. 2005) 
964±54 fBonanos et al. 2006) 


8.2 X 10'^ fSeiEar et al. 2008) 

> 5 X 10^° (Corbelli & Salucci 2000) 


3deg^ 
0.8 deg2 



Table 4. Total amount of data from different survey models (see table [2]). Except for model 7, the surveys are limited 
in Galactic longitude by 20° < 1 < 90°. 



Parameter 


Model 1 


Model 2 


Model 3 


Model 4 


Model 5 


Model 6 


Model 7 


Survey region 


|b|<5° 


|b|<10° 


|b|<5° 


|b|<10° 


lbl<5° 


|b|<10° 


Visible Sky 


Pointing time (s) 


600 


600 


600 


600 


1800 


1800 


120 


Total amount of data (peta-bytes) 


16.7 


33.0 


16.7 


33.0 


50.0 


99.0 


113 



10000 



1000 



100 




10000 



In practice, however, the scattering time can differ from 
this by a factor up to 100 (see Bhat et al. 2004). A smaller 
scattering time leads to the requirement of smaller - and 
thus more - frequency channels. To be on the safe side, 
we apply the highest frequency in the frequency band 
to Eq [5] and divide the resulting scattering time by 100. 
Substituting ^scatt from Eq. [5] in Eq. [7] and taking z^min 
~ 1.115 GHz, leads to a relationship between Ai' and 
DM with a minimum of 42 kHz (at DM = 428). This 
means that having a constant frequency channel width 
of 42 kHz over the frequency band will be sufhcient for 
the pulsar survey. This leads to 9 500 frequency channels 
for i?=400MHz. Figure [5] shows the data rate from a 



Field of the Survey (deg 



Fig. 4. Total amount of data from a FAST pulsar sur- 
vey as a function of the total survey field. The calculation 
assumes a bandwidth of 400 MHz, a FoVbeam of 0.0025 
deg^, 1 polarisation and 2 bits per sample. The three dif- 
ferent lines correspond to observation times of 120, 600 
and 1 800 seconds. The dotted vertical lines indicate from 
left to right the survey fields given by 20° < 1 < 90° and 
|b|<5°, 20° < 1 < 90° and |b|<10° and the entire visible 
sky. 



where iVi 



beams 



is the number of beams in the receiver, t, 



samp 



is the sampling time and Ambits is the number of bits used 
in the digitisation. Aiy can be estimated by demanding 
that the dispersion smearing within the frequency channel 
does not exceed the effective time resolution, given by the 
sampling time and the scattering time: 



Ai/(GHz) 



< 



^samp(Ms) + ^scattC^S 



8.3 X lO^DM 



(7) 



where t'min is the minimum (lowest) frequency in the ob- 
servation frequency band, DM is the dispersi on measure 
and ^ scatt is the scattering time. Following iBhat et al. 
(j2004l ). we relate tscatt to the DM and observing frequency 
as follows: 

logio(iscatt(ms)) = -6.46 + 0.154 logio(DM) -f 

1.07(logioDM)2 -3.861ogio(z^(GHz)). (8) 



CO 

a 




0.01 



FoV (deg'') 



Fig. 5. Amount of data per second from a FAST pulsar 
survey as a function of FoV. The calculation assumes a 
sampling time of lOO/xs, a frequency range of 1.115 to 
1.515 GHz and 2 bits per sample. The dotted vertical lines 
indicate from left to right the FoV of 1, 19 and 100 beams, 
respectively. 



FAST pulsar survey as a function of FoV. Wc have as- 
sumed lOO^s sampling time, a frequency range of 1.115 to 
1.515 GHz and 2 bits per sample. For the initial 19-beam 
receiver, the data rate is just under 0.5 GB per second. 
For 100 beams it is about 2.4 GB per second. 
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The number of operations required to search these data 
for normal pulsars, millisecond pulsars and accelerated 
pulsars in binary systems is approximately 



iVacc X 57Vsampl0g2(A^samp), 



(9) 



where A^dm is the number of DM-trial values and A^acc is 
the number of trial accelerations. Fig. [6] shows the num- 
ber of operations per second required to perform a real 
time analysis of a FAST pulsar survey, assuming 100 trial 
accelerations, a sampling time of 100 /iS, an observation 
time of 600 s and DMmax = lOOOcm^'^pc. A^dm was ob- 
tained numerically by demanding that the total time de- 
lay within the frequency band between two consecutive 
DM-trial values is equal to twice the effective time res- 
olution, given by the sampling time and the scattering 
time of Eq. El This leads to Nt,m = 4 280. For the ini- 
tial 19-beam receiver, the required computation power for 
real time analysis is 9 tera-ops. For 100 beams it becomes 
48 tera-ops. Contributions from de-dispersion, harmonic 
folding and possibly other processes would increase these 
values by a factor less than two. However, it is not unrea- 
sonable to store all the data from a FAST pulsar survey 
and to perform the analysis off-line. This would reduce the 
required computation power significantly. 



Square Kilometre Array comes on-line (§5) FAST will have 
the largest sensitivity of all existing radio telescopes and 
will also be able to see a larger fraction of the sky than 
the Arecibo telescope. Consequently, FAST will be the 
telescope of choice for a number of high-precision tim- 
ing experiments. While the Northern location of FAST 
means that it is unlikely to contribute significantly to the 
timing of, say, the Double Pulsar (PSR J0737-3039A/B, 
Burgay et al. 2003, Lyne et al. 2004), a major contribu- 
tion of FAST will be the participation in the global Pulsar 
Timing Array efforts. 

High-precision timing of a network of pulsars allows de- 
tection of a gravitational wave backgro und by correlatiii g 
the timing residuals of the pulsars fe.g. ljenet et al. 20061 ). 
Current efforts include PPTA (Parkes Pulsar Timing 
Array), EPTA (European Pulsar Timing Array) and 
NANOGrav (North American Nanohertz Observatory for 
Gravitational Waves). In terms of combination of sky 
coverage and sensitivity, FAST will only be matched by 
the European LEAP efforts (Large European Array for 
Pulsars), where the European 100-m class telescopes are 
combined to obtain Arecibo- like sensitivity for pulsar tim- 
ing. This effort, however, has to overcome the challenge of 
coherently combining the signals between telescopes up to 
2 000 kilometres apart. 



1000 




0.001 



0.01 



FoV (deg^ 



Fig. 6. Required computation power for real-time search- 
ing for possibly accelerated fast and slow pulsars in the 
data from a FAST pulsar survey as a function of FoV. 
The calculation assumes a sampling time of 100/is, a fre- 
quency range of 1.115 to 1.515 GHz, 2 bits per sample and 
a DMinax of 1 000. The dotted vertical lines indicate from 
left to right the FoV of 1, 19 and 100 beams, respectively. 



4. Timing 

In addition to being a very powerful search machine, FAST 
will also excel in the high-precision timing of pulsars. For 
a given pulsar, the timing precision roughly scales with 
the S /N of the observed pulse profile (jLorimer fc Kramer 
20051) . and hence with telescope sensitivity. Until the 



4.1. Pulsar Timing Array experiment 



Following e arlier work by Detweilei ( 1979f l and 
Kaspi et al.l (|l994l ). pulsar timing is now routinely 



used to constrain the amplitude of a gravitational 
wave (GW) background. Cur rently the most st ringent 



limits have been placed by Jenet et al. ( 20061 ) from 
observations of seven MSPs with the Parkes 64-m radio 



telescope (cf. ,Manchesteij l2008f) which constrain r^cw: 
the ratio of the energy density of the GW background 
to the closure density to be JIgw^^ < 1-9 x 10~^, 
where h is the dimensionless Hubble constant defined as 
Hq = hx 100 kms^^ Mpc^^. The sensitivity that is needed 
to achieve a GW detection depends on details of the 
Galaxy merger rate, since the GW background generated 
by an ensemble of super-massive black holes distributed 
throughout the Universe is expected to be the strongest 
signal to be detected. The signal originating from such 
sources is expected at a level of flgyvh"^ ."^ ^ ^ 10~^° 
(IJaffe fc Backeij|2003l: IWvithe fc Loebll2003t lEnoki et al 



Simulations suggest that the timing of 20 pulsars 
with a timing precision of 100 ns over five years would 
lead to a sensitivity of ^Iq^^H^ < 6.6 x 10^^^ a nd hence to 



a firs t direct detection of a GW background (j Jenet et al 



2006t ). While this is difficult with current technology. 



FAST will have the sensitivity and sky coverage to make 
a major contribution to achieving this goal. 

Here we estimate the amount of time needed to time 
the brightest MSPs, found by FAST, to high precision. 
For this simulation we choose model 2 from Tabled which 
provides us with 670 MSPs. We take the 50 brightest of 
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these and calculate the amount of time required for a S/N 
of 500, taking into account a conservative profile stabiU- 
sation time of five minutes. This leads to 24 hours to time 
all 50 pulsars once. These numbers do not include accom- 
panying low-frequency observations to allow correction for 
variations in interstellar dispersion. As expected, we found 
that the FoV provided by 19 beams had no impact on the 
timing performance, as there was no instance where two 
or more pulsars could be timed simultaneously. 

5. FAST versus the SKA 

The Square Kilometre Array (SKA) is a planned multi- 
purpose radio telescope with a collecting area approaching 
1 million square metres. Although the exact design is not 
yet determined, it is likely to consist of aperture arrays to 
cover frequencies from 70 MHz up to 500 or even 800 MHz 
and 15-metre dishes for the higher frequencies. About half 
of the elements will be placed within a 5-kilometre core of 
the SKA, the rest will be placed on spiral a rms extend- 



ing se veral thousands of kilometres outwards. ISmits et al 



( 2009[ ) have studied the performance of the SKA concern- 
ing pulsar surveys and timing. Here we will discuss the 
differences between these findings and those of the cur- 
rent paper. 

A pulsar survey with FAST is straightforward com- 
pared to an SKA pulsar survey. With the initial receiver, 
no beam-forming is required. The future PAF does require 
beam-forming, but it does not need the vast computa- 
tional power of several peta-ops required for beam- forming 
the SKA. Another benefit for FAST is that the data do 
not need to be transported over long distances; the data 
processing can be performed close to the dish. Also, the 
data rates and required computation power for a real time 
acceleration search are much lower for the FAST survey. 
This is mostly due to the smaller FoV of FAST, which 
is 10 to 50 times less than the FoV for an SKA pulsar 
survey, but also because these data rates and computa- 
tional power scale with the square of the telescope diam- 
eter. The diameter of FAST is given by the illuminated 
aperture of 300 metres. For the SKA it is determined by 
the largest basel i nes be tween the elements used in the sur- 
vey. Smits et al. ( 2009t ) estimate that the limits in compu- 
tational power will restrict these baselines to 1 kilometre. 
Because of this limitation, the telescope sensitivity of a 
pulsar survey with the SKA is about 2 000m^K~^ which 
is equal to that of FAST. However, the larger FoV of the 
SKA will not only enable a survey of a much larger por- 
tion of the sky, but also allow for a longer dwell time per 
pointing. The actual FoV of the SKA depends on the final 
design, but at the very least it will be equal to the FoV of 
the 15-metre dishes which at 1.4 GHz is about 0.64 deg^. 
For comparison, a reasonable assumption would be that 
an SKA survey has a dwell time of 40 minutes versus 10 
minutes for a FAST survey. This makes the survey sensi- 
tivity of the SKA twice that of FAST. Compared to the 
initial 19-beam receiver of FAST, a pulsar survey with the 
SKA would be about 2.5 times faster, but with the 100- 



beam PAF, the survey speed of FAST would be twice that 
of the SKA. Moreover, because the elements of the SKA 
are fully steerable, the SKA will be able to survey a much 
larger part of the sky and with its latitude of about -30° 
this includes almost the entire Galactic plane. 

Timing pulsars with the SKA can be performed with a 
considerable fraction of the full collecting area, providing 
a sensitivity close to lOOOOm^K"^. Also, the large FoV 
of the SKA allows timing of many pulsars simultaneously 
(but depending on the SKA design, this might not be the 
case when using the SKA for a pulsar timing array). Even 
the future 100-beam PAF of FAST only has a quarter 
of the FoV of 15-metre dishes and should the SKA be 
equipped with PAFs or mid- frequency aperture arrays, its 
FoV increases significantly to 20deg^ or even 250 deg^. 

Given that FAST will be operational well ahead of the 
full SKA, it is likely to be the most powerful telescope 
to perform pulsar science in the next decade and will re- 
main an outstanding instrument for pulsar science, com- 
plementing the SKA once it is operational. 

6. Summary 

The FAST pulsar survey simulation in §2 shows that, de- 
spite the small natural beam-size and limitations in zenith 
angle, FAST will be a formidable instrument for finding 
pulsars. 

A tradeoff between the number of pulsars and total ob- 
servation time can be determined from Table [21 With the 
initial 19-beam receiver about 5 200 previously unknown 
pulsars, including 460 MSPs, can be found in a survey of 
about 230 days of the visible part of the Galactic plane 
with |b| < 5°. With a 100-beam receiver 5 900 previously 
unknown pulsars can be found in the same region in 130 
days, including about 540 previously unknown MSPs. The 
100-beam receiver will also enable an all-sky survey yield- 
ing 4 300 previously unknown pulsars in 300 days. Further, 
we estimate that a 470-hour survey of M31 and M33 would 
yield between 50 and 100 extra-Galactic pulsars. 

A pulsar survey with the 19-beam receiver will produce 
just under 500 MB of data per second. Depending on the 
survey, this leads to a total of 17 or 33 peta-bytes. The 
major part of a real-time acceleration search of these data 
would require 9 tera-ops. After 2014, such data rates and 
processing power should be very feasible. 

FAST will also be able to contribute to the existing ef- 
forts to detect a gravitational wave background by timing 
a large number of MSPs to high precision. Although the 
visible sky of FAST is limited to 58% of the entire sky, 
those pulsars that are visible to FAST can be timed to 
great precision very quickly. A simulation in §4 suggests 
that the 50 brightest MSPs, visible to FAST can be timed 
to a S/N of 500 in just 24 hours. 

Compared to a pulsar survey with the SKA, FAST 
will match the sensitivity of the SKA, but the large FoV 
of the SKA will allow for longer dwell times and make 
the survey faster, at least for the initial 19-beam config- 
uration of FAST. Also, the SKA will be able to survey a 
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much larger part of the sky. However, having a single dish 
rather than many elements spread over a large area signif- 
icantly reduces the required data rates and computational 
power. When timing pulsars for a pulsar timing array, the 
SKA will have much more sensitivity than FAST, since 
the SKA can use almost the full collecting area in this 
case. Also, the larger FoV of the SKA will allow timing of 
many pulsars simultaneously, although this will probably 
not benefit the timing of the limited number of pulsars 
from the timing array. Given that FAST will be opera- 
tional well ahead of the full SKA, it will provide the best 
prospects for pulsar science in the next decade. 
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